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Abstract 

The baroclinic instability problem is considered in the framework of Laplacian tidal theory. The 
Hilbert space of the quasigeostrophic vorticity budget is spanned by spheroidal functions. The fluid 
is linearly stable against quasigeostrophic disturbances. As the essential source of irregular ocean- 
atmosphere motions, baroclinic instability is ruled out by tidal theory. The midlatitude /3-plane 
budget of vorticity fluxes is inconsistent with basic laws of motion on the rotating spherical surface. 
Realistic numerical simulations of global wave dynamics and dynamical circulation instabilities 
require a covariant account of fluid motion on the spherical planet. 

Pacs numbers: 47.32.-y, 47.35. +i, 92.10.Hm, 92.60.Dj 
I. Introduction 

Major difficulties in numerical climate simulations result from the fact that the potential energy 
in the ocean-atmosphere system is closely associated with large-scale features of its density field 
while a considerable fraction of its kinetic energy resides on fairly small eddy scales [|l| . State-of-the- 
art models for weather prediction and climate simulation capture large-scale features of the global 
circulation with some degree of realism if they are determined by large-scale features of topography 
and external forcing. However, simulation of the dynamical control of density by transfer processes 
between the small scales of kinetic energy and the large potential energy scales is less satisfactory |2j . 
The key to these energy- and vorticity-fluxes is fluid instability. Of the numerous instabilities in the 
climate system none is considered as fundamental as baroclinic instability ||. This concept refers to 
the growth of weakly divergent Rossby waves in a stably stratified and vertically sheared fluid on the 
rotating spherical surface. The contemporary understanding of cyclogenesis, predictability-limits and 
the transition to chaos, turbulence and stochasticity in the climate system revolves essentially around 
this process. Hence, it also provides the paradigm for the design and interpretation of irregular fluid 
motion in numerical circulation models, ranging from resolved scales down to the parameterization of 
subscale processes. 

Baroclinic instability theory invokes four sets of approximations. First: shallow water theory 
with constant mean layer-thicknesses and -velocities Q. Second: quasigeostrophy. For tidal theory, 
Longuet-Higgins j|| has shown that this approximate closure of the vorticity budget yields a meaningful 
and elegant Rossby wave filter of Laplace's tidal equation. In baroclinic fluids this approach governs the 
thermal wind. Third: the midlatitude /3-plane. To avoid formal difficulties with spherical coordinates, 
a tangential cartesian plane is pinned to the sphere at some extratropical latitude in analogy to 
Kelvin's well-established equatorial /3-plane. Fourth: a number of scaling assumptions which neglect 
fluid velocities relative to planet rotation and stratification relative to barotropicity wherever it appears 
uncritical. Quantitatively, these assumptions seem generally to be justified ||. 

The paper at hand considers the baroclinic instability problem of Laplace's tidal theory, indepen- 
dent of the midlatitude /3-plane and scaling assumptions. Tidal theory differs in two aspects from 
current circulation theories and models: it is Newtonian (namely: covariant) and accounts consis- 
tently for the globe's sphericity Q. The relationship of its analytical structure and the theory of Heun 
functions |7j is increasingly understood. While a transformation of the tidal equation into Heun's 
equation is not known, general approximations and exact special cases can be expressed in terms of 
Heun functions. These concepts provide a sound foundation for Rayleigh stability analysis in the 
framework of global wave-circulation theory. On the basis of covariant fluid dynamics on the rotating 
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spherical surface it will be shown that Rossby wave growth in baroclinic fluids plays a far lesser role 
for irregular ocean-atmosphere motions than currently thought. 

The geometrically and dynamically consistent way to avoid formal difficulties with spherical co- 
ordinates is the use of index notation and covariant differentiation. Here, indices m, n, ... = 1,2 run 
over longitude A and latitude ip while braced indices (£) = 1,2 refer to the top and bottom layer of 
the fluid and are not subject to the summation convention. Covariant differentiation will be denoted 
by a semicolon. For details of the notation see ||. In this formalism, equations in curvilinear coordi- 
nates look widely similar to corresponding equations in cartesian coordinates with geometrical details 
consistently absorbed into core symbols and indices. Thus, formulas emphasize the physical structure 
of the problem. 

II. Tidal Equations 

The problem is considered in terms of the bishallow water equations on the rotating spherical 
surface || . Rayleigh stability theory requires the linearization of such equations around the considered 
basic state. For both layers, the mean layer mass per unit area and hence the concentration are here 
assumed to be constant 

R = Rn\ + Rr 2 ) = const, r = R^/R = const. 
Thus, the effective pressure P(R, r) and the interfacial potential u(R, r) 

p = o 7(2) (1 + 5r 2 )R 2 , fi = 7 (2) S Rr, 
are also constant where 7^ = g/pm, g the gravitational acceleration and 

S = (P(2) -P(i))/P(i) > 

the positive definite stratification parameter. The spherical generalization of a constant mean barycen- 
tric velocity V n and a constant mean vertical shear W n are given by 

V n = a 2 U (cos 2 <p, 0), W n = a 2 W(cos 2 ip, 0) 

with Earth's radius a and constant angular velocities 

U = rU (1) + (1 - r)U {2) , U x = (1 - r)U (1) + rU {2) , W = U (1) - U {2) . 

The barycentric and baroclinic potential vorticities of this circulation , Zq and Zi, are obtained as 

RZ = F = 2(Q + U )sinip, RZ X = F X = 2(tt + U^sinip 

while S = Wsinip. This circulation is driven by an external, meridionally varying surface pressure 
and the value of the vertical shear W is determined by f2, the stratification parameter 5 and the 
equator-to-pole gradient of the surface pressure Q. For small-amplitude perturbations (rn,Tj, jm^n) 
of the state vector: layer mass, concentration, barycentric and baroclinic mass flux, linearization of 
the bishallow water equations around this circulation leads to the tidal problem || 



d m + j n ; n =-RW n d n n (1) 

Rd lV + i n ; n = -r 12 W n d n m (2) 

dojn + e mn F j m + d nP * = -W m 8 m i n - 2e mn Si m (3) 

dii n + e mn F x i m + d ni i* = -r 12 (W m d m j n + 2e mn Sj m ) (4) 
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where r\ 2 = r(l — r) and d /i = &t + Uq/i&x- The linearized pressure emerges as 

P* = (<9j? -P) r m + (<9 r P)j? rj = c 2 pim + c 2 p 2 Rr] 
while one finds for the linearized interfacial potential 

/i* = ri 2 R(dR /j,) r m + ri 2 R(d r h)r rj = c^im + c 2 fi 2 Rrj. 

Here, c 2 = 7(2) -R arid 

Pi = 1 + 5r 2 , nx = r/i2 = r i2 p2 = 5rrx2- 

These coefficients satisfy 

c 2 (pi + M2) = c 2 (l + p 2 ) = Co + c 2 , c 4 (pi^ 2 - P2M1) = c 4 /i2 = CqC 2 

where the intrinsic barycentric and baroclinic phase speeds are given by 



c o/i = ^c 2 (! + Sr±y/(1- 5rf + 45r 2 ). 

Taking the curl of (3) and (4) one arrives at the perturbation vorticity budgets 

R 2 d z + fd a F = -W n d n R 2 ( - 2i a 8 a S (5) 

R 2 d!( + i a d a F l = -r 12 (W n d n R 2 z + If 3 a S) (6) 
with barycentric perturbation vorticity 

R 2 z = e an j n ; a -F m-2SRrj (7) 

and baroclinic perturbation vorticity 

R 2 ( = e an i n]a -F 1 Rr ] -2r 12 Sm. (8) 

Equations (1) through (8) pose the Rayleigh stability problem for a generic stably stratified and 
vertically sheared fluid on the rotating spherical surface. In spherical bishallow water theory, constant 
mean layer-thicknesses exclude a mean flow with available potential energy. The tidal problem for a 
circulation with finite available potential energy is discussed in the Appendix. 

III. Quasigeostrophic Stability Analysis 

In the barotropic 1-layer limit, equations (1) through (4) reduce to Laplace's standard tidal equa- 
tions. In special cases, exact analytical solutions of the tidal equation are known in terms of confluent 
Heun functions, namely spheroidal functions @j §]■ In the entire wave number space, approximate 
analytical solutions can be expressed in terms of spheroidal functions and the asymptotics of tidal 
functions coincide with the asymptotic behaviour of prolate spheroidal functions: the Margules regime 
of globally defined Legendre polynomials at small Lamb parameters ||] and the Matsuno regime of 



Hermite polynomials on the equatorial /3-plane at large Lamb parameters [10]. Longuet-Higgins' 
quasigeostrophic Rossby wave filter Q retains these functional characteristics. In the 1-layer limit, 
the perturbation vorticity budget becomes 

R 2 d t z + fd a F = 

with perturbation vorticity 

R 2 z = e an j n ; a -Fm. (9) 
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As in the strictly nondivergent case, quasigeostrophy assumes that the perturbation mass flux is 
sufficiently represented by a stream function A 

A — f pim, A 
Jn — c nm u 

while the mass perturbation in (9) is not supposed to vanish, thus giving rise to the notion of weakly 
divergent perturbations. A closed expression for (9) is now obtained by invoking the geostrophic 
approximation 

c 2 m = -FA 
leading for the vorticity budget to the equation 

(A - aV - Mt)A = 

where y = sirup, a = 2a(Q + U)/c the Lamb parameter, M the zonal wave number, v = a(u — UM)/c 
the Doppler-shifted frequency and r = ajv. This is the prolate spheroidal wave equation. The 
dispersion relation for quasigeostrophic Rossby waves becomes 

v = -aM/e(N,M;a) 

with prolate spheroidal eigenvalue e(N, M; a). Comparison of this expression with numerical solutions 
of the complete tidal equation (fig.l) demonstrates that with the exception (of the gravity branch) 
of the Yanai wave (mode number N=0) quasigeostrophy provides a satisfactory approximation to all 
Rossby modes of the tidal problem. Also, this shows that quasigeostrophy is by no means a regional, 
e.g. extratropical approximation. Rather, it is globally valid and includes the Margules regime as well 
as the Matsuno regime. Physically, the equatorial /3-plane approximation of the prolate spheroidal 
equation accounts for wave trapping in the Yoshida guide. A similar wave guide in midlatitudes 
does not exist and a midlatitude /3-plane does not appear in the systematic approximation theory of 
spheroidal functions ||. 




Figure 1: Tidal cigenfrcqucncies (solid lines) and quasigeostrophic approximation (dotted lines). Frequencies 
larger than \M\ correspond to gravity modes. Negative zonal wave numbers M indicate westward propagation. 
Dashed-dotted line: v = a. Tidal frequencies were calculated by [11] and spheroidal eigenvalues computed with 
NAG-Lib routine F02GJE. 
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The functional structure of the tidal problem uniquely determines the physical interpretation of the 
so-called "beta-effect" and the spectrum of the tidal wave operator. The Coriolis term of tidal theory 
represents the meridional shear of a mean zonal flow with uniform angular velocity f2 + U. Doppler- 
shifts only appear with respect to U since the observer corotates with f2. While the corotating observer 
does not see a frequency shift with respect to f2, the corresponding meridional shear does not vanish 
on this transformation. Physically, the "beta-effect" of tidal theory refers to such meridional shear 
and differs profoundly from topography on the f-plane. 

The spectrum of the tidal wave operator represents free waves in an elastic medium. Such a medium 
has two types of excitations: longitudinal P (primary, pressure or sound) waves and transversal S 
(secondary or shear) waves. Long gravity waves of covariant shallow water theory are represented as 
longitudinal sound waves with restoring vertical buoyancy forces appearing as an effective "compress- 
ibility" of the strictly 2-dimensional system. Correspondingly, shallow baroclinic gravity waves are 
represented as second sound in a strictly 2-dimensional bifluid. Rossby waves, on the other hand, obey 
the dynamics of (radially polarized) shear waves 12]. The frequencies of S-waves are always lower than 



P-wave frequencies: all Rossby frequencies lie below gravity frequencies. Rossby waves are transver- 
sal and mean flow shear exerts restoring stresses tangential to wave fronts. Low-frequency Rossby 
waves are essentially divergence- free and well represented by Margules' approximation. High Rossby 
frequencies are limited by meridional trapping at low latitudes which induces a weak divergence: weak 
divergence is characteristic of equatorial Rossby waves. Tropically trapped, weakly divergent Rossby 
waves are well approximated by Matsuno's theory. Longuet-Higgens' quasigeostrophy unifies both 
approaches. With two Lame coefficients, elastic wave theory is inherently a two-parameter problem. 
This is also true for generic tidal theory: while a "compressibility" controls gravity wave dynam- 
ics, the Lamb parameter represents the mean meridional shear that governs Rossby wave dynamics. 
Quasigeostrophy is the shear wave filter of tidal theory. 

In application to the present bishallow problem, quasigeostrophy represents the barycentric and 
baroclinic mass flux perturbations in terms of stream functions 

a — e pflri \ ■ _ , r) m ')h 

Jn — c nm u <"n — c nm u Y 

and determines mass- and concentration-perturbations in (7) and (8) from the thermal wind relation 
of (3) and (4) 

(?pi m + c 2 p 2 Rrj = — FqA — 2Sip 
c 2 [ii m + c 2 [i 2 Rf] = —Fill) — 2r\ 2 SA. 
Solving for m and Rrj and inserting the result into (7) and (8) the perturbation vorticities become 

R 2 z = -(A - h y 2 )A - hy 2 ^ 

R 2 (=-(A-h iy 2 )iP-r 12 hy 2 A 

with 

h = (M2«(2) + r i2a? 2 )// i 2, h = (rp 2 a 2 2) + «i)/a«2, h = (P2«( 2 ) - a 12 ai)/n2 

where am = 2a(f2 + Uf£\)/c, a\ 2 = am — a( 2 ) and a\ = 2a(fl + Ui)/c. With these expressions the 
vorticity budgets (5) and (6) assume the form 

(u; A - H 00 y 2 - f M)A = (WM(A + 2) - H m y 2 )^ (10) 

(wiA - H u y 2 - /iM)^ = r 12 (WM(A + 2) - H 10 y 2 )A. (11) 
Here, o; /i = u ~ Uq/\M and / /i = 2(0 + Uq/i) while 

H 00 = u h + r 12 hWM, H m = u; h + h x WM 
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and 

#10 = uih + h WM, Hu = ujihi + r 12 hWM. 

Equations (10) and (11) are a system of coupled spheroidal equations and spheroidal functions form 
a complete set of eigensolutions. Eliminating Aip from these equations yields 

qrl) = (A - h y 2 - Mt )A 

with r w = 2(17 + U^))/(lo - U^M) and 

q = MT 12 -hy 2 , t 12 = r (1) - t (2) , r = rr (1) + (1 - r)r (2 ), n = (1 - r)r (1) + rr (2 ). 

Substituting this expression for the baroclinic stream function into (10) results in a single fourth-order 
equation 

(A - h iy 2 - Mn)-(A - /i y 2 - Mr )i = r 12 <^4 (12) 

q 

for the barycentric stream function. Utilizing now the spheroidal property: A.A = ((3 2 y 2 — e)A with 

e + Mr - (/3 2 - h )y 2 = Xq 
for constant X, equation (12) requires the simultaneous validity of the two quadratic equations 

t 12 (X 2 + (1 - 2r)X - r 12 ) =0 

X 2 h + (h - hi)X - r l2 h = 0. 

These two equations express the major difference between the cartesian and the spherical stability 
problem. In cartesian geometry, the dispersion relation determines admissible eigenfrequencies. On 
the rotating spherical surface, background inhomogenities due to the planet's sphericity, coordinate- 
dependent Coriolis forces and the mean circulation also determine admissible Lamb parameters. 
Rayleigh theory of spatially inhomogeneous systems accounts for wave trapping. The compatibil- 
ity of both quadratic equations is determined by the relation 

fy) - hi = (1 - 2r)h - «(i)a( 2 )//x 2 (13) 

which follows from the definition of ho, hi and h. Given (13), essentially two classes of solutions 
of (12) exist: either, the stratified background is free of vertical shear or Rossby waves are strictly 
nondivergent. In the first case: W = and t\ 2 = 0. Hence: To = t\ = r, while q = —hy 2 and (12) 
becomes 

(A - f3 2 y 2 - Mr)l(A - [3 2 y 2 - Mr) A = 

where the Lamb parameters are obtained as 

(3 = 2a(n + U)/c , f3i=2a(n + U)/ci 

with intrinsic barycentric and baroclinic phase speeds c /i- The dispersion relation is 

(uo e(N, M; ft) + ftM)(^ e(N, M; ft) + ftM) = 

with prolate spheroidal eigenvalue e(N,M;f3) and vq/i = a(uj — UM)/co/i- In the shear-free case, 
Rossby waves are weakly divergent, propagate as barycentric and baroclinic modes and eigenfrequen- 
cies are real. In the second case: c — > oo and ho = hi = h = 0, while q = Mn 2 = const so that (12) 
reduces to 

(A-Mr ( i))(A-Mr (2) )A = 0. 
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In the strictly divergence-free case, spheroidal functions degenerate into Legendre polynomials and 
the dispersion relation becomes 

(e + Mr (1) )(e + Mr (2) ) =0 

with e(N, M) = N(N + 1) + (2N + 1)|M| + M 2 . The mean vertical shear traps nondivergent Rossby 
waves in individual layers and eigenfrequencies are real. This remains true if the mean circulation 



exhibits available potential energy [13|. 



There are two more solutions if in one of the layers a«) = 0, i.e. Um = — 0. For au\ = 0, Rossby 
waves exist only in the complementary layer £ ^ j and from (12) one finds the dispersion relation 

e(N,M;f3 {e) ) = -Mr {i) 

with Lamb parameter 

(3f f) =Aa 2 (n + U (e) ) 2 /g'H {e) 

where Hm = Rt(\/pn\ is the mean layer-thickness and g' = g(pm ~ P(i))// 9 (2) the reduced gravity 
accelaration. Both of these solutions are weakly divergent and stable. The condition U(j\ = —Q 
implies a westward flow circulating the globe in one day. Hence, the layer is at rest in a nonrotating, 
inertial system and such conditions do not admit Rossby wave propagation (the mean meridional shear 
is absent). In the rotating system, corresponding large-scale velocities are of a magnitude that is not 
met on this planet. For practical purposes these solutions are hence of little significance. 

The same dispersion relations follow if A and ip in (10) and (11) are replaced with spheroidal 
functions and the resulting system of 3 algebraic equations is solved for the unknown Lamb parameters 
and eigenfrequencies (see Appendix). As a further alternative, the layer representation of spherical 



linearized bishallow water [ 13 1 may be chosen as starting point for the stability analysis rather than 
the modal representation (1) through (4). It is readily seen that the dispersion relation from such an 
approach coincides with the results obtained above (see Appendix). The baroclinic instability problem 
of tidal theory does not assume the form of a spherical Taylor-Couette flow and all eigenfrequencies are 
real. Unlike baroclinic gravity waves, Rossby waves do not feed on the energy of a stably stratified and 
vertically sheared mean flow. Thus, the flow is linearly stable against quasigeostrophic disturbances. 

For the system (1) through (4) isopycnals coincide with equipotential surfaces. This type of 
configuration is generally considered in baroclinic instability theory: isopycnals are assumed to be 
"flat" and a slope-parameter does not enter the problem ||. On the other hand, observers and 
modellers are typically concerned with sloped isopycnals and their erosion by baroclinic instability 
P]. The tidal equations for bishallow water with sloping isopycnals are well known and their stability 



against nondivergent Rossby waves has been demonstrated [13]. The stability of such a system against 
quasigeostrophic Rossby waves is shown in the Appendix. 

The physical interpretation of these results is best considered in comparison to Kelvin-Helmholtz 
instability. For this instability, highly divergent baroclinic gravity waves continuously sample both 
layers of a stably stratified and vertically sheared fluid. The dynamics of baroclinic gravity waves are 
controlled by the competition of (stabilizing) stratification and (destabilizing) vertical shear. If the 
vertical shear becomes too large, baroclinic gravity waves grow. None of these mechanisms plays a role 
in tidal Rossby wave dynamics. Rossby waves are inseparably linked to a definite value of the mean 
meridional shear determined by O + Urg\ . Hence, they propagate in individual layers and a distinction 
of barycentric and baroclinic Rossby waves is meaningless (unless the vertical shear vanishes). Weak 
divergencies are associated with meridional trapping and do not accommodate the exploration of 
adjacent layers. None of the mean flow features of a stably stratified, vertically sheared flow competes 
with the restoring stresses. Thus, a transfer of energy or vorticity between this type of circulation and 
Rossby waves is excluded in the framework of tidal theory. 

Although mean flow available potential energy may alter the system's wave guide geography rad- 
ically, it does not change its stability properties. For very low-frequency Rossby waves, isopycnal 
slopes modify the effective Lamb parameter (see Appendix). This parameter may now assume real or 
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imaginary values. In the event of imaginary Lamb parameters, spheroidal wave operators change from 



prolate to oblate || and Rossby waves may be meridionally trapped in a polar wave guide |0|, 14, 15]. 
Although the Lamb parameter of oblate Rossby waves becomes imaginary, their frequencies remain 
real j|] and Rossby wave amplitudes do not commence growing. In general, a stably stratified and 
vertically sheared fluid with or without available potential energy lacks the faculties of energy- and 
vorticity-transfer to Rossby waves. This statement may require modification for certain initial condi- 
tions or the extremely steep isopycnals associated with outcropping. Independent of the role, these 
and other special cases may take, Laplace's tidal theory does not support the ubiquity, baroclinic 
instability is currently thought to have for irregular ocean-atmosphere motions. 

IV. Discussion 

The results of the previous section are in clear contrast to baroclinic instability theory on the 
midlatitude /3-plane. To identify the source of this discrepancy evaluate the midlatitude /3-plane 
approximation of (12). In this sense, the latitude y is fixed at some (extratropical) value while 
Gq = hoy*, G\ = h\yl and G 2 = hy 2 are considered as constant parameters. Adopting furthermore a 
cartesian Laplace operator and trigonometric eigenfunctions with K 2 = a 2 {k 2 +k2) and M = ak\ cost/?* 
one obtains from (12) the approximate dispersion relation 

(K 2 +G\+ Mt x )(K 2 + G% + Mr ) = r 12 q 2 . 

A little algebra readily shows that this expression is equivalent to 

(K 2 + F 2 1} + Mt {1) )(K 2 + F 2 2) + Mr (2) ) = sF^Ffo 

with s = P(i)/p(2) and F 2 ^ = P 2 ^y 2 - This equation is quadratic in frequency. Unlike the results 
of the previous section, its roots do not assume an easily interpreted form. This is indicative of 
the difficulties of the primarily geometrical midlatitude /3-plane approximation to accommodate the 
system's physical structure. Nevertheless, the reality of these roots is obvious and the midlatitude 
/5-plane approximation to the final wave equation (12) comes qualitatively to the same result as the 
spherical analysis with respect to stability. Hence, the trigonometric approximation of wave functions 
in itself - though unsatisfactory - is uncritical. The source of discrepancies is therefore a fundamentally 
different account of vorticity fluxes by covariant shallow waters and baroclinic instability theory. 

In the covariant case, the vorticity budget (12) is uniquely determined by the equations of motion 
(1) through (4) and the quasigestrophic thermal wind approximation. The validity of quasigeostrophy 
for the barotropic fluid is well demonstrated by fig.l and the particular form of the same argument for 
the baroclinic system is again an unambiguous consequence of the equations of motion. These equa- 
tions are the direct and unique result of the application of covariance requirements to the formulation 
of spherical shallow water dynamics. Thus, the vorticity fluxes of covariant shallow water theory are 
in essence the expression of basic geometrical and physical consistency conditions for the hydrostatic 
flow of the stably stratified and vertically sheared fluid on the rotating spherical surface. 



Standard baroclinic instability theory, on the other hand, takes the Primitive Equations [16] as a 
starting point and invokes the midlatitude /3-plane approximation to derive the vorticity budget || . It 
has been shown that the Primitive Equations do not pose a covariant dynamical problem and involve 
ambiguous mass and momentum fluxes as a consequence of the violation of Newton's first law §. 
For the vorticity budget it now becomes important that the midlatitude /3-plane is much more poorly 
defined than the equatorial /3-plane. Effectively, it takes the role of a geometric closure assumption 
which introduces ill-defined vorticity fluxes. In the stability analysis this inconsistency resurfaces as 
spurious Rossby wave growth, i.e. baroclinic instability. 

Vorticity fluxes are crucial in maintaining and changing the density field of the ocean-atmosphere 
system and a consistent representation of vorticity dynamics is indispensable for realistic numerical 
simulations of the global circulation. At this time, a large number of numerical circulation models is 
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based on the Primitive Equations. Moreover, most models are not formulated in terms of spherical 
coordinates but utilize a multi-/3-plane approach to approximate the globe's sphericity: Laplace oper- 
ators, for instance, are coded as a sum of second order derivatives, ignoring first order contributions 
from nontrivial Christoffel symbols. Widely independent of spatio-temporal resolution and the quality 
of subscale parametrizations, such models cannot expect to simulate the large-scale circulation with 
geometric-dynamic integrity and realism. 

Theoretically as well as numerically, large-scale ocean-atmosphere dynamics require a covariant 
dynamical framework including the acknowledgement of the spherical geometry of the planet's surface. 
The linear stability of the baroclinic fluid against quasigeostrophic disturbances calls for a covariant 
reanalysis of the observational evidence on dynamical instabilities in the ocean-atmosphere system. 



Appendix 

The effect of available potential energy of the circulation is considered. With constant surface 
pressure, the mean layer mass per unit area is now assumed to vary with latitude according to 

R m = R%\l-b {t) y 2 ) 

I 

water equations, its slope parameters 



where R E is its equatorial and R p its polar value. As a geostrophic solution of the nonlinear bishallow 



>{t) 

satisfy 

6 (1) = a 2 (2n + U {1) + U {2) )W/2 1 R { z ) « cfnW/jR^. 
b {2) = a 2 [(2Q + U (2) ){5U {2) -W)- U (1) W]/2 7 Rf « a 2 n(5U (2) - W)hR%\ 

with 7 = 7(2)5. The associated mean vertical shear W = — Uq) is in thermal wind balance and 
the Coriolis parameter in layer £ is given by 

F(t) = R(t)Z{i) = 2(n + U^)) sirup. 
With dm = dt + U^d\, the layer representation of this tidal problem assumes the form [13] 

d(e)m {e) +jfyy, n = 

d {i) v n ] + + d n P{£) = 

where the perturbation pressures are given by 

P(l) = 7(1)^(1) + 7(2)™(2), P(2) = 7(2) ("1(1) + m {2)) = 7(2) "I- 

The curl of the momentum budgets yields the vorticity budgets 

R(e)d(e)Z(£) + jfydaZy) = 

for the perturbation vorticities 

R(£)Z(£) = e an v^; a -Zmm^y 
The quasigeostrophic approximation with 

A i] = e nm d m A (e) 

and the thermal wind relation 

7m(i) = Z(2)A(2) - z (i) A (i), S7m( 2 ) = sZ (1 )A (1 ) - ^(2)^(2) 
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leads to the coupled wave equations 

7 [(#(i)A - R a (l) d a - Ffo/^d^ - e an R 2 {1) d a Z {1) d n }A (1) = -R {l) Z {2) F {l) d {l) A {2) 

7 [(i?( 2 )A - Rf 2) d a - Ff 2) / S1 )d {2) - e an R 2 (2) d a Z (2) d n ]A {2) = -R {2) Z (1) F {2) d {2) A {1) . 

For < 1 6(f) | <C 1 these equations are a coupled system of spheroidal equations. Physically, this 
condition on the slope parameters excludes outcropping of isopycnals. With 

0fa = 4a 2 (n + U^f/g'Hf, h {i) = (3 2 {e) + 2b {i) Mr {i) 

the quasigeostrophic vorticity budgets assume the form 

(A - h {1) y 2 - Mr (1) )i (1) = -sk(3 {1) (5 {2) y 2 A {2) 

k(A - h {2) y 2 - Mr (2) )A {2) = -(3 {1) (3 {2) y 2 ' A {1) 

where k = \] H$ /H$. Note that will become negative for very low Rossby wave frequencies 
at negative b^M. Using now AA^ = (f3 2 y 2 — ejAw, these equations reduce to a system of coupled 
algebraic equations for the amplitudes of the stream functions. Nontrivial solutions exist for vanishing 
coefficients of the polynomial X + Yy 2 + Zy 4 . For the present problem, these coefficients have the 
form 

X = (e + Mr (1) )(e + Mr (2) ) =0 

Y = (e + Mr {1) )(/3 2 - h {2) ) + (e + Mr (2) )(/3 2 - h {l) ) = 

Z = (3 4 - (h {1) + h {2) )(3 2 + h (1) h {2) - sp 2 {1) (3 2 2) = 

The first of these equations states the reality of eigenfrequencies for stream functions which remain 
regular at the poles. Also, it demonstrates that Rossby waves propagate in layers unless the vertical 
shear vanishes. The second equation selects the Lamb parameter for the respective mode and the 
third equation determines the Lamb parameters. For = the solutions of the main text emerge 
while the case c — > oo and b^ ^ has been discussed in fl3| . 
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